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By using observables that only depend on charged particles (tracks), one can efficiently suppress 
pile-up contamination at the LHC. Such measurements are not infrared safe in perturbation theory, 
so any calculation of track-based observables must account for hadronization effects. We develop a 
formalism to perform these calculations in QCD, by matching partonic cross sections onto new non- 
perturbative objects called track functions which absorb infrared divergences. The track function 
Ti{x) describes the energy fraction a; of a hard parton i which is converted into charged hadrons. We 
give a field-theoretic definition of the track function and derive its renormalization group evolution, 
which is in excellent agreement with the Pythia parton shower. We then perform a next-to- 
leading order calculation of the total energy fraction of charged particles in e^e~ — >■ hadrons. To 
demonstrate the implications of our framework for the LHC, we match the Pythia parton shower 
onto a set of track functions to describe the track mass distribution in Higgs plus one jet events. 
We also show how to reduce smearing due to hadronization fluctuations by measuring dimensionless 
track-based ratios. 



Jets are collimated sprays of particles that arise from 
the fragmentation of energetic quarks and gluons. Nearly 
every measurement at the Large Hadron Collider (LHC) 
involves jets in some way, either directly as probes of 
physics in and beyond the Standard Model, or indirectly 
as a source of backgrounds and systematic uncertainties. 
In order to predict jet-based observables using quantum 
chromodynamics (QCD), one typically performs infrared- 
and collinear-safe (IRC safe) jet measurements which in- 
volve only the kinematics of the jet constituents [T]. In 
particular, IRC safe jet measurements do not distinguish 
between charged and neutral particles, despite the fact 
that, for example, charged pious (tt^) are measured us- 
ing both tracking and calorimetry whereas neutral pions 
(tt") are measured using calorimetry alone. 

In this letter, we develop the theoretical formalism to 
calculate track-based observables, which depend on the 
kinematics of charged particles alone but not on their in- 
dividual properties or multiplicities. The experimental 
motivation for track-based measurements is that track- 
ing detectors offer better pointing and angular resolution 
than calorimetry. By only using tracks, one can sub- 
stantially mitigate the effects of pileup (multiple colli- 
sion events in a single bunch crossing) which is becom- 
ing more relevant as the LHC achieves higher luminosity 
(see e.g. [2H5] for alternative approaches). In addition, 
tracks can aid in jet substructure studies where the an- 
gular energy distribution in the jet discriminates between 
different jet types [H [8] . While we focus on charged par- 
ticles, this formalism applies to any (otherwise) IRC safe 
measurement performed only on a subset of particles. 

To describe track-based observables in QCD, we in- 
troduce the track function Ti{x,fi). A parton (quark or 
gluon) labelled by i with four-momentum hadronizcs 
into charged particles with total four-momentum = 
xp^ + O(Aqcd). The distribution in the energy fraction 



< a; < 1 is the track function and is by definition nor- 
malized 

/ dxT,{x,^Ji)^\. (1) 
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The track function is similar to a fragmentation function 
(FF) or a parton distribution function (PDF) in that 
it is a fundamentally non-perturbative object that ab- 
sorbs infrared (IR) divergences in partonic calculations. 
Like FFs and PDFs, the track function has a well-defined 
dependence on the renormalization group (RG) scale fi 
through a DGLAP-type evolution [5HI31 , though the spe- 
cific evolution is more reminiscent of the jet charge dis- 
tribution [m [TS] . The track function description is valid 
for factorizable observables where partons are well sepa- 
rated and can be considered to hadronize independently. 
Corrections to this approximation are captured by power 
corrections, which are beyond the scope of this letter. 

Consider the cross section for an IRC safe observable 
e measured using partons 

where we drop possible convolutions with PDFs to keep 
the notation simple. Here, Hn denotes A^-body phase 
space, da N / dllN is the corresponding partonic cross sec- 
tion, and e{{pi}) implements the measurement on the 
partonic four-momenta . Since e is an IRC safe observ- 
able, the KLN theorem [16l|T7] guarantees a cancellation 
of IR divergences between real and virtual diagrams. The 
cross section for the same observable measured using only 
tracks is 
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FIG. 1: LO (dotted) and NLO (solid) track functions ex- 
tracted in Pythia from the fraction of the jet energy carried 
by charged particles. 



where Ti{xi) is the track function for parton i. In the 
absence of track functions, this cross section would con- 
tain collinear divergences due to a mismatch between real 
and virtual diagrams. The track functions absorb these 
IR divergences, and the partonic cross section is cor- 
respondingly modified with respect to cr at. We will show 
below how to perform this matching explicitly for the 
example of e^e^ — > qqg. The universality of collinear di- 
vergences in QCD p!8H20] guarantees their cancellation 
to all orders in the strong coupling constant as- 

At leading order (LO) in as , the cross section depends 
on a single partonic multiplicity N and there are no IR 
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divergences implying a 
simply a finite distribution which can be obtained di- 
rectly from the energy fraction of charged particles in a 
jet initiated by a parton i. In Fig. [T] we show the track 
functions obtained from pure quark and gluon jet samples 
produced by Pythia 8.150 [52] and clustered using 
the anti-fcr algorithm [23^ in Fast Jet 2.4.4 [24|. (To ex- 
tract the track function at next-to- leading order (NLO) 
we use Eq. (11); the jet radius R is correlated with the 



RG scale /i.) As expected, the up- and down-quark track 
functions are very similar, with a peak at x = 0.6. This 
means that on average 60% of the energy of the initial 
quark is contained in charged hadrons, in agreement with 
a recent CMS study [53]. The small difference between 
up and down is due to strangeness, since us mesons are 
charged whereas ds mesons are neutral. Because gluons 
have a larger color factor than quarks, they yield a higher 
track multiplicity, and the corresponding track functions 
are narrower, as expected from the central limit theorem. 

Formally, the (bare) track function is defined in QCD 
in a fashion analogous to the unpolarized FF (cf. ^51127) ). 
Expressed in terms of light-cone components. 



Tqix) = J dy+ d y±e 
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FIG. 2: The evolution of the NLO gluon (top) and d-quark 
(bottom) track functions compared to Pythia. Starting from 
II — 100 GeV (shown in Fig. [TJ, we evolve using Eq. ^ down 
to ^ = 10 GeV and up to /i = 1000 GeV. 



where tp is the quark field, C (N) denote charged (neu- 
tral) hadrons, and is the large momentum compo- 
nent of all charged particles. Whereas the FF describes 
the energy fraction carried by an individual hadron, the 
track function describes the energy fraction carried by 
all charged particles. As for the FF, the definition can 
be made gauge invariant with eikonal Wilson lines. The 
gluon track function is defined analogously [21] • 
At NLO, the bare track functions bare 

in pure di- 



mensional regularization with d = 4 — 2e are 
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which arise from collinear splittings, controlled by the 
timelike Altarelli-Parisi splitting functions Pi^jk{x) 1^. 
Renormalizing the ultraviolet divergences in Eq. ^ m 
MS leads to the evolution equation for the track function 

^J.-^Ti{x,^l) ^l^y^ [dzdxidx2^^^^Pi^jk{z) (6) 

d/i 2^ J TT 

3,k 

X Tj{xi, fl)Tk{x2, ^J-) S[x~zxi-{l-z)x2]- 

Like for a PDF, the track function can be extracted at 
one scale and RG evolved to another scale, and the evo- 
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lution preserves the normalization in Eq. ([T]). Unlike a 
PDF, Eq. ^ involves a convolution of two track func- 
tions at NLO (and more convolutions at higher orders 
corresponding to multiple branchings), so it is numeri- 
cally more involved to perform the /i-evolution. At lead- 
ing logarithmic (LL) order, the RG evolution in Eq. ^ 
is equivalent to a parton shower, and Fig. [2] demonstrates 
excellent agreement between our numerical evolution and 
the parton shower in Pythia. 

For a calculation at NLO, both the partonic cross sec- 
tion and the track functions have IR divergences which 
cancel in Eq. ^ . To demonstrate this in a simple exam- 
ple, consider the process e^e~ — > hadrons at a center- 
of-mass energy Q where one measures the total energy 
fraction w of charged particles. At NLO the partonic 



process is e~^e~ — ^ qqg, whose kinematics are described 
by the energy fractions yi and y2 of the quark and anti- 
quark. Applying Eq. ^ we find 

= [dyidy2 , [ dxidx2dx3Tq{xi)Tq{x2)Tg{x3) 

dw J dyidy2 J 

X S{w - [yixi + y2X2 + (2 - ?/i - y2)x3]/2). (7) 
At the LO partonic level 
da(0) 
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where is the total Born cross section. At NLO, the 
cross section can be expressed using plus-functions as 
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where Cp — 4/3. The l/eip-divergences in a^^'' are can- 
celled by the ones in T(i) from Eq. and the finite 
remainder defines a^^'' in MS. (For different IR regula- 
tors, the track function may also contribute finite terms 
to the matching.) By performing this kind of matching 
calculation, one can determine ctjv for any process to any 
order in a^. For the case of Eq. ^ which involves a sin- 
gle scale Q, we choose /i ~ Q to minimize the logarithms 
in a. In Fig. [3] we show the LO and NLO distributions 
for the energy fraction w, using the track functions ex- 
tracted from Pythia by means of Eq. (111. There is 



good convergence from LO to NLO and our fixed-order 
calculation agrees well with Pythia. 
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FIG. 3: Normalized distribution of the energy fraction w of 
charged particles in e^e~ at Q = 91 GeV, calculated at LO 
(green) and NLO (orange), compared with Pythia (blue). 
The uncertainty bands are obtained by varying fi between 
Q/2 and 2Q, and do not include track function uncertainties. 



Ultimately, we are interested in applying the track 
function formalism to jet-based measurements at the 
LHC, which typically involve multiple scales. For nar- 
row well-separated jets, though, the energy fraction x of 
the charged particles only depends on the jet scale 



/ij ~ ptR , 



(10) 



where pt is the transverse momentum of the jet and R de- 
notes its azimuthal-rapidity size as defined by a jet algo- 
rithm (anti-ZcT in this letter). Indeed, by varying R (and 
trying different jet algorithms) while keeping fij — pxR 
fixed, we find nearly identical x distributions in Pythia. 
At LO, the normalized distribution in x defines T^°\x) 
itself, shown in Fig. [T] At NLO, the cross section has the 
same form as for the jet charge distribution |15) 
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and we can invert Eq. (Ill to determine the NLO track 
functions in Fig. [T] The jet functions Ji and coeffi- 

I and 



cients Jij were determined at NLO in Ref. 
Refs. [m 130] , respectively. The matching coefficients Jij 
first appeared for fragmentation inside an identified jet in 
Refs. [SUini], which discuss their physical interpretation. 

As a track-based measurement relevant for the LHC, 
consider the track-only jet mass spectrum in pp to Higgs 
plus one jet. This example is more complicated than the 
e^e~ example above since it involves several scales: the 
hard scale set by the transverse momentum of the jet pj.. 
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FIG. 4: Track mass distribution in pp H+]et obtained 
from the Pythia parton shower matched onto either track 
functions or the Lund string modeL 
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FIG. 5: Comparison between track and calorimeter measure- 
ments of ratio of jet mass and jet pr- Here the cuts are on 
and f)j from the tracks in the jet. 



the jet mass scale mj, and the scale associated with soft 
radiation TOj/p^. In a resummed jet mass calculation, 
one would need to treat each emission in the exponenti- 
ated soft function as hadronizing independently, i.e. con- 
voluted with its own track function |28j . At LL order, 
however, we can use the fact that the parton shower al- 
ready describes the numerous parton emissions that build 
up the jet mass distribution. Thus, a correct and instruc- 
tive use of our formalism at LL is to run the perturbative 
Pythia parton shower down to a low scale /x ~ Aqcd 
and match each final state parton onto a track function 
(also evolved to ^ ~ Aqcd)- 

In Fig. [4] we show the track mass toj spectrum us- 
ing the E-scheme [SS] for the anti-/cT jet algorithm with 
R— 1.0. We impose realistic cuts on the jet rapidity rjj 
and pj. at the calorimeter level (using all the particles) , 
to select the same events for both curves. There is excel- 
lent agreement between the track function and the Lund 
hadronization model [M] in Pythia, which is a non- 
trivial check considering that the track functions were 
originally extracted from Pythia using jet energies (and 
not jet mass) at /Lt = 100 GeV. Since the track function 
does not include correlations between the hadronization 
of different partons, the agreement in Fig. |4] shows that 
independent fragmentation is an excellent approximation 
for these observables. The small difference between the 
two distributions in the peak region is due to nonpertur- 
bative power corrections, which to first approximation 
can be described by a shift in w?j (see e.g. [551 155H57] ) . 

Despite the advantages of using track-based observ- 
ables for pile-up suppression, it should be noted that the 
track functions have a rather large width (see Fig. [ij. 
Fluctuations in the charged energy fraction produce an 
effective energy smearing for a track-based measurement 
relative to a calorimetric one. To partially address this 
issue, one can focus on observables for which the track 
function would only have a modest effect, such as for di- 
mensionless ratios of observables. A particularly useful 
example are iV-subjettiness ratios |38j . which are relevant 
for jet substructure studies. Though a calculation of such 



ratios is beyond the scope of this letter, we plot in Fig. [5] 
the ratio of the jet mass to jet pr in Pythia, measured 
using either tracks alone or all particles (calorimeter) . As 
expected, the smearing effect is reduced since hadroniza- 
tion fluctuations are correlated between the numerator 
and denominator. 

Given the potential experimental benefits from track- 
only measurements at high luminosities, we expect that 
track functions will offer an important theoretical han- 
dle for future precision jet studies at the LHC. Beyond 
the tests of our formalism against Pythia performed 
here, we stress that it is possible to systematically im- 
prove the accuracy of track-based predictions. At fixed 
order in as, one can calculate higher-order corrections 
for the (process-dependent) matching onto track func- 
tions. More ambitiously, one could match track func- 
tions onto automated NLO calculations, perhaps using a 
more convenient IR regulator (such as dipole subtraction 
|39j ) rather than the MS scheme used here. To improve 
the accuracy of resummation, one needs to determine 
higher-order RG evolution of the track functions, and 
precision track-based studies would also require power 
corrections. Ultimately, one would want to follow the ex- 
ample of PDFs and extract track functions for the LHC 
from global fits to data. 
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